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Abstract. Measurements of the surface brightness distri- 
bution and of the velocity dispersion profile have been so 
far used to infer the inner dynamics of globular clusters. 
We show that those observations do not trace back the 
dark matter potentially concealed in these systems in the 
form of low-mass compact objects. We have built Michic 
models of globular clusters which contain both massive 
and low-mass stars. An analytic expression for the stel- 
lar mass densities has been explicitely derived in terms of 
the usual error function and Dawson's integral. While the 
heavy population is kept fixed, the abundance of the light 
species of our models is varied. When stellar velocities are 
anisotropic, both the surface brightness and the velocity 
dispersion profiles of the cluster become insensitive to the 
abundance of low-mass stars. This suggests that the ac- 
tual stellar mass function of many globular clusters is still 
to be discovered. 



1. Introduction 

Among the different kinds of astrophysical star conglom- 
erates, globular clusters (GC) have a very peculiar status. 
First, as we go down the mass-scale of observed struc- 
tures, from super clusters through galaxies down to dwarf 
spheroidals, they are the first objects whose dynamics can 
be understood without resorting to the presence of dark 
matter. Second, and this may be related to the first point, 
they are the only star clusters dense enough for a sta- 
tistical equilibrium to be reached in a Hubble timescale. 
The dark matter content of globular clusters may be an 
important clue to the understanding of the formation of 
our galaxy. However, the observational situation is not 
very clear yet as regards stellar counts and the related 
direct detection of a hidden population of light and faint 
stars. Since the claim by Richer and Fahlman that the 
mass function steeply rises towards low masses in some 



GCs (Fahlman et al. 1989, Richer et al. 1991), the Hubble 
Space Telescope (HST) has been used to measure the mass 
function in several globular clusters, and the opposite con- 
clusion was reached (Demarchi and Paresce 1995, Elson et 
al. 1995, Paresce, De Marchi and Romaniello 1995). Part 
of the discrepancy is due to the conversion scheme used 
to translate the measured luminosity functions into mass 
functions. There is actually no consensus on the mass-to- 
luminosity relation of low-mass objects. In a recent inves- 
tigation, Santiago et al. (1996) found a steadily increasing 
mass function for u> Cen, using HST observations. This 
suggests that the dark matter content may depend on the 
cluster. It seems therefore important to fully understand 
what we can also learn from other types of measurement, 
such as in particular the brightness and the velocity dis- 
persion profiles. 

There is a widespread belief that globular clusters can- 
not contain large amounts of dark matter, because of the 
Virial theorem. The latter relates the velocity dispersion a 
at the center of the cluster to the total mass M t and to the 
half-mass radius ?*h of the system. For a one-component 
globular cluster, that relation may be expressed as (Spitzer 
1987) : 

.GM t (i) 
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However, if a dark component is also present, in the form 
of low-mass objects for example, the virial theorem must 
be modified to 



(2) 
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where the index i refers to the different stellar species, and 
where E p is the potential energy of the cluster. In previous 
papers (Taillet et al. 1995 and 1996) focusing on isotropic 
globular clusters, we showed that in the presence of dark 
matter, relation ([!]) still holds for the visible component, 
and that most of the visible properties of the cluster are 
approximately the same as for a one-component system. In 
particular, the brightness profile is nearly unchanged even 
if dark matter is abundant. This is due to the phenomenon 
of mass segregation. Because of gravitational interactions, 
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stars tend to share their kinetic energy. Thus, the light ob- 
jects gain velocity from the heavier stars, which sink to the 
center of the cluster, while the light objects have a much 
more diffuse distribution. The heavier stars then form a 
self-gravitating system in the low-density background of 
the light objects, and their properties are not significantly 
altered. It was concluded that even a dark-matter domi- 
nated cluster looks very much the same as a dark-matter 
free system, and that the point of the dark matter con- 
tent could not be made by the mere observation of the 
brightness profile. 

It was also noted that the presence of dark matter 
would be betrayed by the velocity dispersion profile of 
bright stars which flattens when dark matter is present. 
However, those profiles need to be accurately measured, 
which is not yet the case. The results of our analysis 
were obtained for King models where the stellar distri- 
bution is isotropic in velocity space. Their validity must 
be reassessed when the velocity distribution is no longer 
isotropic. Anisotropy is likely to be present at some level 
because stars mainly interact in the cluster core. At least 
part of those orbiting at the outskirts of the system gained 
the required energy by interacting at the center. Then, 
orbits should be predominantly radial at large radii. In 
that respect, note that oj Cen for which the HST observa- 
tions provide a steadily rising mass function may also have 
an anisotropic velocity distribution (Meylan 1987), even 
though an isotropic distribution is not excluded (Merrit 
et al. 1996). Such a cluster may be modeled by a Michie 
distribution in phase space, where an angular-momentum 
exponential cut-off is applied to a King distribution (see 
section 2). 

It is well known that when anistropy is present, the ve- 
locity dispersion profile is less fiat. Thus we suspect that 
the flattening of the velocity dispersion profile that oc- 
curs in the presence of dark matter will be partially can- 
celled by the effect of anisotropy. A given velocity dis- 
persion profile could therefore be interpreted either with 
a one-component isotropic GC or with a two-component 
anisotropic GC. However, anisotropy may also affect a pri- 
ori the brightness profile. The main goal of this paper 
is therefore to make these statements more quantitative. 
We will investigate how anisotropy affects the detectabil- 
ity of dark matter in globular clusters. In section 2, we 
present a new analytic form for the equation that relates 
the gravitational potential to the stellar mass densities. 
The Michie distribution is expressed with the usual error 
function and Dawson's integral. Section 3 is devoted to our 
set of two-component Michie models whose heavy compo- 
nent is kept fixed while the abundance in low-mass stars is 
varied. Various velocity anisotropies are explored. Finally, 
the surface brightness and the velocity dispersion profiles 
of these models are discussed in section 4 and conclusions 
are drawn. 




Fig. 1. The unit vector e s is aligned along the line of sight 
HM that crosses the globular cluster. At the pericluster H, 
the distance to the center O is minimal. At any given point 
M on the line of sight, with position labelled by the angle 
77, a convenient frame may be defined as presented here. 



2. The Michie distribution functions 

A Michie model of a globular cluster is spherically sym- 
metric as regards its distribution of stars in space. How- 
ever, its velocity ellipsoid tends to be elongated towards 
the center of the cluster for radii larger than some critical 
value r a . The phase space distribution of a stellar species 
with mass m and whose one-dimensional velocity disper- 
sion is a, may be expressed as 

f(r,v) = ke- L2 / 2r > 2 {e^> 2 - l} , (3) 

where L = OM A v is the orbital momentum of the stars 
located at point M and whose velocities are v. That orbital 
momentum is defined with respect to the cluster center 
O. The quantity £ is related to the stellar energy per unit 
mass E = $(r) + v 2 /2 through 

S = $ t - E , (4) 
where $ t = &(rt) is the gravitational potential at the 
tidal boundary r t of the cluster. Whenever £ is negative, 
the distribution function / vanishes. Inside the radius of 
anisotropy r a , the distribution of stellar velocities is spher- 
ical. Beyond r OJ it straightens in the direction of the clus- 
ter center. At remote distances, trajectories are almost ra- 
dial. This is a major difference with King models for which 
velocities are spherically distributed. The mass density is 
now given by a double integral. The orbital momentum 
L = r v sin 8 depends actually on both the magnitude v of 
the velocity and on the angle 9 between the latter and the 
radial direction OM. Defining the variables x — v/\/2a 
and y = cos 9 yields a mass density 

p{r) = 47rfc m (2cr 2 ) 3/2 H(u,a) (5) 
proportional to the double integral 

rl r\fu 

H(u,a) = dy dxx 2 h(x,y) , (6) 
Jo Jo 
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where the function h(x, y) is defined as 
h( X ,y) = e -«V(l-y 2 ) [ e (u-x 2 ) _ 1 



(7) 



for the dispersion. In this article, the dispersion velocity 
5 ra d along the line of sight will be defined as 



The parameter u = (x — ip) is the difference between the 
reduced potential \ — ^t/c 2 at the tidal boundary r t of 
the system and its counterpart ip = <J>(r)/cr 2 at distance 
r from the cluster center. The ratio r/r a is denoted by a. 
The scaling factor between the mass density p(r) and the 
integral @ may also be inferred from the relation 

P(r) = H{u,a) 

where p(0) denotes the stellar mass density at the cluster 
center. 

So far, expression (|^) was estimated numerically through 
a direct integration. We have derived here an analytic de- 
velopment of Tl(u, a) in terms of the functions 



«rad 



2^ N 



1/2 



(14) 



d(x) 
and 
e(x) 



dt 



dt 



(9) 



(10) 



The first integral is called the Dawson's function and may 
be easily computed thanks to a convenient approximation 
due to Rybicki (Rybicki 1989 and Numerical Recipes). 
The function e{x) is related to the error function 



e(x) = 



erf(x) 



(11) 



The distribution function /(r, v) can be averaged over the 
angle 9 to yield a King distribution up to a correction 
factor 

£/cj 2 _ x \ \d{ax) 



f(r,v) = k \ e 



(12) 



The parameter ax stands for the ratio rv / \2r a <j . Inside 
the radius of anisotropy, ax is small compared to unity 
and the correction factor d(ax) / ax of relation ( |T^ ) reduces 
to 1. A King distribution is therefore recovered when the 
radius r is smaller than r a . After an integration on the 
magnitude v of the velocity, relation (||) may be expressed 
as 



H(u, a) 



1 



2(1 + a 2 ) 



d(a-\fu) — ay 



(13) 



Several approximations to that expression are given in the 
Appendix depending on whether yfu and a^fu are large 
or not with respect to unity. 

The dispersion velocity profile of globular clusters may 
be defined in various ways. What astronomers actually 
measure is the velocity of some stellar species along the 
line of sight. Suppose that a specific pixel contains N stars 
with same spectral type. The radial velocity of the i th ob- 
ject of the sample is denoted by Uj. A variety of statistical 
averages can be performed, each yielding a different value 



In our two component models, we will be interested in the 
velocity dispersion profile of the heavy and bright stellar 
population for which spectral measurements can be per- 
formed. 

Let us consider now a line of sight with direction set 
by the unit vector e s as shown in Fig. |l|. At the pericluster 
H, the distance to the center O is minimal. For any given 
point M along that line of sight, a convenient frame may 
be defined with its unit vector e z pointing outwards in 
the same direction as the vector OM. The unit vector e x 
is perpendicular to e z and is in the plane defined by the 
three points O, H and M. Finally, the last basis vector e y 
is defined as the product e z A e x . The position of point 
M may be traced by the angle r\ between the vectors OH 
and OM. The projection v s of the velocity along the line 
of sight is then given by 

e s ■ v = v x cos rj + v z sin 77 . (15) 
Its square v 2 may be averaged locally to yield 
(v 2 s ) = (v 2 x )cos 2 V + (« 2 ) sin 2 77 . (16) 
Because the velocity distribution is axisymmetric around 
the direction connecting point M to the cluster center O, 
the contribution of the product v x v z to the local average 
of the line of sight velocity v 2 vanishes. The latter may be 
expressed in terms of the functions I(u, a) and J{u : a) 

p(M) ((^y^=4irkm (2a 2 f 2 {l + J sin 2 77} (17) 

where p(M) denotes the mass density of the stellar species 
with mass m under consideration at the point M. The 
functions X and J are respectively defined by the integrals 

l(u,a) = dy dxx A (1-27 2 ) h(x,y) , (18) 
Jo Jo 

and 

J{u,a) = I dy I dxx 4 (3y 2 - l) h(x,y) . (19) 
Jo Jo 

The functions TL and X are related by 

l(u,a) = "0 , (20) 

so that X may be expressed as 

3 

(l + a 2 ) X(u, a) = H + — g {d (ay/u) — a\/u\ 



(21) 



The calculation of the function J is slightly more involved 
and yields 

(l + a 2 ) 2 J(u,a) = ^{e(V^)-V^} 



(2 2 \ ^ 

a + u + ua ) 



3/2 



/C(it, a) {d (a\fu) — a^/u] , (22) 
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where 



K,{u, a) = 



2a 2 



U + a 2 ) (5 + 2ua 2 



4a 5 



(23) 



Notice that v 2 needs eventually to be averaged along the 
entire line of sight. The end result is therefore the radial 
average 



are assumed to have a standard solar L/M ratio. Low 
mass objects are too faint to shine and their luminosity is 
negligible. 

The structure of the cluster is determined by the Poisson 
equation where the local mass density is given by 



p(r) 



Pel 



H(u, a) 



Introducing the dimcnsionlcss radius z — r/a where the 
typical scale length a is defined as 



p{M) (v 2 ) ds 



7i(m2u/mi, a) 



(25) 



■rail 



(24) 



p(M) ds 

to be compared with the observational definition (|l~4|). 



leads to the differential equation 



(26) 



3. The models 



z 2 dz 



, du 
dz 



TL{u, a) 

n( x ,o) 



Pc2 
Pel 



Ti.(rri2u/mi, a) 
Ti(m 2 x/mi,0) 



To test the sensitivity of observable quantities to the amount Once that equation is solved and the structure of the clus- 



of dark matter, we built two-component Michie models of 
globular clusters, with both a heavy and a light stellar 
population. One of the components encapsulates the vis- 
ible, solar mass stars while the other component is dark 
matter in the form of low-mass objects with mass rri2 = 
0.1 Mq. The indexes 1 and 2 respectively refer to the 
bright and dark populations. We built 25 globular clusters, 
with anisotropy radii z a = 1000, 100, 50, 20 and 10, and 
dark matter content M2/M1 — 0, 1, 3, 6 and 10. The pa- 
rameter z a stands for the ratio of the anisotropy radius r a 
to the typical scale length a yet to be defined. Models with 
z a = 1000 may be considered as King models. In that case, 
the anisotropy radius is so large that it contains most of 
the cluster, hence a spherical velocity distribution almost 
everywhere. A two-component Michie model is completely 
specified by six parameters, namely the anisotropy radius 
r a , the velocity dispersions a\ and o 2 of the two species, 
the normalization factors k\ and &2 (or alternatively the 
central densities p c \ and p c2 ), and finally the depth of the 
gravitational potential well \ = $ t / ' o\ . These parameters 
are entirely determined by the following conditions : 

— (i) The velocity dispersion for the luminous stars is 
measured. We will adopt the typical value o\ = 7 km/s. 

— (ii) The two species are in thermal equilibrium, impos- 
ing the relation m\<j\ = TO2<j| which sets 02 equal to 
VIOxcti. 



ter is determined, the velocity dispersion and the surface 
brightness profiles are computed. The former is evaluated 
with the help of relation (|24|). The latter is just the surface 
mass density of heavy stars 



Si = / Pl (m)ds , (27) 
up to the factor (Lq/Mq)/4it. 



4. Discussion and conclusions 

Profiles of the surface mass density Si and of the veloc- 
ity dispersion v v:i d for the heavy component of our models 
are presented in Fig. |^ and |^ respectively. Low mass stars 
do not shine in the V band and therefore do not con- 
tribute to the surface luminosity of the cluster. Figure 2 
shows that whatever the anisotropy radius z a , the bright- 
ness profiles are almost not affected by the dark matter 
content of the cluster. At constant z a , the various curves 
are superimposed while the light to heavy component mass 
ratio M2/M1 increases from to 10. This result was al- 
ready obtained in the case of King models (Taillet et al. 
1995 and 1996) and its validity is now extended to Michie 
clusters. For z a = 10 and 20, the surface luminosity drops 
like r~ 5 / 2 as expected for a strongly anisotropic system. In 
Fig. |^, the velocity dispersion profiles are presented. For z a 
(iii) The visible central density is set to p cl = 8000M Q /pcL 1QQ ^ velocities are almost spherica ii y distributed and 
This leads to models having almost the same central the curveg are different when the dar k matter abundance 

M2/M1 varies. However, the velocity dispersion along the 
line of sight becomes insensitive to the dark matter con- 
tent of the cluster when the anisotropy of the velocity 
distribution increases. For z a smaller than 20, the various 
profiles start to be fairly similar and become hard to dis- 
entangle observationally. Beyond the radius of anisotropy 
r a , the stellar orbital velocity vg — v sin 6 is dominated 
on average by the Michie exponential cut-off on the or- 
bital momentum. The orbital velocity is in that case of 
order vg ~ ar a /r. It is no longer sensitive to the gravi- 
tational potential. The velocity dispersion v ra( j measured 



brightness which, as a matter of fact, is a well mea- 
sured quantity. 

— (iv) The anisotropy radius has been varied as explained 
above. The smaller r a , the stronger the anisotropy. 

— (v) The total visible luminosity is also a well deter- 
mined quantity. Here, it has been set equal to a typical 
value Li = 3 x 1O 5 L . Throughout our set of models, 
p C 2 and x have been chosen to yield that total lumi- 
nosity, as well as the above mentioned ratios M2/M1. 

Remember that Mi and M2 respectively stand for the 
total mass in heavy and in low mass stars. Visible stars 



Frangois Roueff et al.: Velocity dispersion : a closer look 



5 




Fig. 2. In each of the four plots, the surface mass density 
Ei of the heavy population of our models is presented 
as a function of the radius r. The mass ratio M2/M1 is 
set equal to the five different values 0, 1, 3, 6 and 10. 
The radius of anisotropy z a decreases from 1000, where a 
King model is recovered with a spherical distribution of 
velocity, down to 10 for which the anisotropy is strong and 
the stellar orbits almost radial. In each plot, the curves 
are hard to disentangle. The surface brightness profile is 
insensitive to the dark matter content M2 of the cluster. 

along the line of sight, to which the orbital velocity con- 
tributes most, is therefore blind to the mass content of 
the cluster. As z a decreases, the effect is more and more 
pronounced. When it becomes dominant everywhere out- 
side the core radius of the cluster, the gravitational po- 
tential and therefore the dark matter of the system have 
little effect on u ra d, hence the degeneracy of the profiles 
in Fig. ||for z a = 10. Note that as velocities become radi- 
ally distributed, light stars need more space to contribute 
the same dark mass as a result of the conservation of the 
phase space volume. At fixed Ma jM\ ratio, low mass stars 
extend further away as the velocity anisotropy increases. 
However, there is still a substantial amount of them in- 
side the visible part of the cluster. When z a = 10 for in- 
stance, the dark mass inside the inner 60 pc is respectively 
4.7 x 10 4 , 1.3 x 10 5 , 2.2 x 10 5 and 2.9 x 10 5 M for a light 
to heavy stellar mass ratio M2/M1 of 1, 3, 6 and 10. 

The dashed curve that appears in three of the plots of 
Fig. U is the velocity dispersion profile for a King model 
with no low mass stars. That case is degenerate with the 
entire set of Michie clusters with z a = 20, whatever their 
dark matter content. The same velocity dispersion profile 
may therefore be interpreted either with a one-component 



Fig. 3. Same as in figure 2 with the velocity dispersion 
profiles. For z a = 1000, velocities are almost spherically 
distributed and the curves are different when the dark 
matter abundance M2/M1 varies. However, the velocity 
dispersion along the line of sight becomes insensitive to the 
dark matter content of the system when the anisotropy of 
the velocity distribution increases. The dashed curve that 
appears in three of the plots is the velocity dispersion pro- 
file for a King model with no low mass stars. That case 
is degenerate with the entire set of Michie clusters with 
z a = 20. The same velocity dispersion profile may there- 
fore be interpreted either with a one-component isotropic 
cluster or, alternatively, with an anisotropic system con- 
taining dark matter. 

isotropic cluster or, alternatively, with an anisotropic sys- 
tem containing dark matter. If so, the degeneracy may 
be lifted by the distribution of surface luminosity. Note 
however that the brightness profiles become different far 
from the center, in a region where the surface mass density 
Si has already dropped by four orders of magnitude with 
respect to its central value. Figures 2 and 3 illustrate actu- 
ally the difficulty to estimate the amount of low mass stars 
potentially concealed in globular clusters from their veloc- 
ity dispersion profiles and from the distributions of their 
surface luminosities. Both measurements have been used 
so far to determine the mass content of GCs. We there- 
fore conclude that the same set of observations may be 
interpreted by quite different models. This suggests that 
the actual stellar mass function of many globular clusters 
needs to be reinvestigated. 

Dark matter in globular clusters should be traced by 
the three dimensional velocity distribution of bright stars. 
We have just shown that the projection along the line of 
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sight alone does not lead to a unique answer. Another pos- 
sibility is the search for tidal tails. If present, dark matter 
dominates the gravity of the outskirts of clusters and pre- 
vents heavy stars from evaporating (Moore 1996). Tidal 
tails exist whenever dark matter is not present to play that 
inhibition effect. As mentioned by Taillet et al. (1995 and 
1996), the presence of low mass stars should also lead to 
an infrared halo surrounding the bright central part. Fi- 
nally, if the cluster lies against a rich stellar background, 
low mass stars should also induce a few gravitational mi- 
crolensing events. 
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or ( AS ) should be used in expression (|13|) . If both variables 
are small at the same time, the expansion becomes 



u 5/2 

H{u,a) = — AiA 2 {l + A 3 u(l-a 2 ) 

+ A 3 A 4 u 2 (l - a 2 + a 4 ) 

+ (l - a 2 + a 4 - a 6 ) 

+ A 3 A A A b Asv 4 (l - a 2 + a 4 - a 6 + a 8 ) 

+ ...} 

In the same limit, the function X may be expanded as 



(A3) 



I[u, a) 



,7/2 



-A 1 A 2 A 3 {l + A 4 u(l~2a 2 ) 



+ A 4 A 5 u 2 (1 - 2a 2 + 3a 4 ) 

+ A A A 5 A 6 u 3 (l - 2a 2 + 3a 4 - 4a 6 ) (A4) 
+ AiA 5 AQA 7 u 4 (l - 2a 2 + 3a 4 - 4a 6 + 5a 8 ) 
+ ...} 

whereas J is a sum of three series expansions 
,2 , u 9 ' 2 



(l + a 2 ) z J(u,a) 
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■a 2 A 3 Ai {2J a (u) 



- (5 + 7a 2 )a 2 J b (ua 2 ) 

+ 9(l + a 2 )a 2 J c (ua 2 )} . (A5) 
The functions J a , J b and J c are respectively defined by 

J a (x) = 1 + A 5 x + A 5 A 6 x 2 + A 5 A 6 A 7 x 3 + ... (A6) 

J b (x) = 1 - A 5 x + A 5 A 6 x 2 - A 5 A 6 A 7 x 3 + ... (A7) 
while 

J c {x) = 1 - A 4 xJ b {x) . (A8) 



Appendix A: Useful expansions 

The functions d(x) and e(x) may be expanded as power 
series of the variable x when the latter is small. The Daw- 
son's function is given by 

d(x) = x - A ± x 3 + AxA 2 x 5 + ... 

+ (-l) n A 1 A 2 ...A n x 2n+1 + ... , (Al) 

whereas 

e(x) = x + Aix 3 + A x A 2 x h + ... 

+ A^.AnX 2714 - 1 + ... . (A2) 
The constants A n stand for the ratios 2/(2n + 1). When- 



ever cty/u or y/u are small compared to 1, expansions (Al) 



